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Abstract 

Due to strong interactions between superfluid particles, many particles are not in the condensate 

o 

CN \ near absolute zero. In nonequilibrium situations produced under external stirring, they contribute 

to the flow via their induced average macroscopic motion. We present a modification of the Bose 

hj ; 

If} . condensate model that considers this induced internal motion. Specifically, we show that the 

<N . 

velocity field of the internal motion acts as a gauge field with which the condensate interacts, 
' thereby coupling the applied vorticity to the superflow. The resulting Ginzburg-Landau-like model 

o 



predicts the flow phenomena observed near absolute zero, in particular, transition to turbulence. 
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I. INTRODUCTION 



A proper description of superfluid hydrodynamics must account for transition to turbulent 
flow. Turbulent excitations, in general, are generated by external 'stirring' mechanisms that 
produce vorticity. They occur in the superfluid component in the form of atomically thin 
vortex filaments of quantized strength, a complex tangle of which forms the turbulent state. 

At T £=: 0, the fundamental microscopic theory of superfluids, based on many-body quan- 
tum mechanics, is difficult to use to study flow phenomena. The macroscopic model of a 
weakly interacting Bose condensate provides a simpler approach. It captures some essential 
physics of superfluids, in particular, the structure of the vortex lines. Continuum mechanical 
equations of superfluidity are usually built on the assumption that the condensate model 



gives the exact description at zero temperature [1]. The main drawback of this model, 
however, is that superfluid is a liquid and not a dilute gas. Because of strong particle in- 
teractions, a large number of superfluid particles are not in the condensate at T « Q 
and these contribute to the flow in nonequilibrium situations. We refer to these excited 
particles as the zero-point cloud, to be distinguished from the thermal cloud (normal fluid) 
at finite temperatures. Our basic aim is to rectify this drawback by a modification of the 
condensate model that considers the contribution of the zero-point cloud. We show that in 
nonequilibrium situations produced under external stirring, an average macroscopic motion 
of the cloud is induced; whose velocity field provides a new degree of freedom for the liquid. 
This velocity acts as a gauge field with which the condensate particles interact, entailing 
a gauge coupling between the interpenetrating flows of the condensate and the cloud. The 
applied vorticity, thus, couples to the superflow in a gauge-field manner by means of in- 
ducing macroscopic internal vorticity via the cloud. The resulting model coincides with 
the Ginzburg-Landau (GL) model of superconductivity [3j in the superconductive phase, 
and explains the peculiar similarity between the latter and superfluid hydrodynamics, in 
particular transition to turbulent flow. 

II. THE EQUILIBRIUM CONDENSATE MODEL 

In the standard condensate model, the large number of condensate bosons per unit vol- 
ume, no, gives rise to a collective irrotational velocity field, v = V0, for the condensate, 
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which corresponds to the superfTuid velocity. The macroscopic wave function of the con- 
densate is given by ip = -Juq exp(im<fi/ti), where m is the boson mass. For a condensate 
of bosons interacting via a delta-function repulsive potential of strength g, the equilibrium 



energy is given by the Gross-Pitaevskii energy functional 

d 6 x 
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where /i > is the chemical potential near absolute zero. The field equation that follows from 
minimization of this functional _f the time-independent Gross-Pitaevskii equation) possesses 
vortex line solutions (see e.g. [4j); the core being an atomically thin hole surrounded by a 
macroscopic region of irrotational flow. Thus, under external stirring, formation of vortex 
lines becomes favorable in appropriate flow geometries and the superfTuid can lower its total 
energy by creating inhomogeneity in the flow, concentrating vorticity along highly localized 
regions. 

The Goldstone mode associated with the spontaneous breaking of the global gauge sym- 
metry of (00) corresponds to the long-wavelength perturbations of the velocity potential. It 
plays the main part in the excitation spectrum near absolute zero. The liquid in the (sound) 
wave moves irrotationally with velocity given by the gradient of the perturbation field. Ev- 
idently, such a long-wavelength mode cannot exist in a disordered state, which renders the 
model inadequate for the description of turbulent flow. However, as we shall see, with the 
appearance of additional interactions in nonequilibrium situations, the long-wavelength ex- 
citations become costly and the Goldstone mode disappears. This is where gauge potential 
enters, as we explain in the sequel. As an inspiration, we first remark below the intriguing 
correspondence between superfluidity and superconductivity. 



A. Superfluidity versus superconductivity 

As for classical fluids, a superfluid undergoes transition from laminar to turbulent flow 
with or without intermediate mixed (vortex) flows, depending on the flow geometry. Ex- 
amples of these are the cylindrical vessel (or Couette) flow and pipe flow, respectively (see 
js] and the references therein). In the superfluid vessel flow, as long as the angular velocity 
of the rotating vessel is below a first critical value the fluid remains stationary. The 
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intermediate mixed flow forms when the angular velocity is increased above fii, with the 
superfluid breaking into an ordered array of vortex lines that are aligned along the rotation 
axis. Above a second critical angular velocity ^2 (> ^1), the vortex flow develops into a 
turbulent flow in which quantum vortices invade the whole fluid by forming a chaotic tangle. 
In pipe flow, however, turbulent flow sets in when the superfluid velocity exceeds a critical 
value and there is no intermediate mixed flow. Except for the formation of vortex lines, the 
condensate model can explain none of these hydrodynamic behaviors. 

The intriguing correspondence of the situation here with that in superconductivity is 
inspiring for our gauge-field description. External stirring agents weaken the ability of a 
flow to remain laminar, just as external currents weaken the ability of a metal/alloy to 
be superconductive. In both examples above, expulsion of the applied vorticity from the 
liquid interior corresponds to the Meissner effect, which occurs below a critical value of the 
applied magnetic field. The effect breaks down when the applied vorticity is too large. The 
ensuing flow can be divided into two types according to how this breakdown occurs, which 
is determined by the flow geometry. The flow types correspond to types I superconductivity 
and II. The array of vortices in the mixed superfluid vessel flow corresponds to the Abrikosov 
lattice of magnetic flux lines. The latter forms the mixed state of type II superconductors 
when the magnetic field lies between two critical values. If the field increases even more, 
the flux lines will invade the whole sample and superconductivity will be lost. In type I 
superconductivity, which corresponds to the pipe flow, this happens when the external field 
exceeds a critical value and there is no intermediate mixed state. 



III. GAUGE INTERACTION AND THE GINZBURG-LANDAU MODEL 

In an imperfect Bose gas, the fraction of particles in the condensate is not unity. This 
fraction is even smaller in the case of a Bose liquid with strong particle interactions [2[. That 
is, the densities of the condensate (mn ) and of the superfluid are not the same. However, 
the superfluid velocity is the same as the condensate velocity, v , which is determined by the 
condensate particles alone. The particles in the cloud do not act cooperatively and, hence, 
their atomic momenta cannot contribute to the collective irrotational flow. Nonetheless, 
they do give rise to effective atomic currents, the current density of which is a rapidly 
fluctuating quantity. This entails velocity fields that vary appreciably on the atomic scale of 
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dimensions, the average of which over a macroscopic volume only being pertinent. Physically, 
in view of the randomness of these velocity fields, one expects zero average and hence no net 
macroscopic velocity contribution. But what if the cloud particles are subject to external 
stirring? To still suppose that the induced motions are completely random and not correlated 
with the applied vorticity is counterintuitive. In nonequilibrium situation established under 
the influence of external stirring, an average internal motion must be induced and, in fact, 
this is how the external vorticity couples to the superflow. The condensate model does not 
take this coupling into account, which renders it inadequate for a proper description of flow 
phenomena, in particular, turbulence. Let us, therefore, envisage that the externally induced 
microscopic motion of the large number of cloud particles, each with its atomic momentum, 
gives rise to an average macroscopic internal velocity field V. (This is analogous to the 
induced internal magnetization due to effective atomic currents that result from an applied 
magnetic field, and does not entail cooperative behavior, of course.) The velocity fields 
V and v are regarded as independent variables so that the liquid is endowed with a new 
degree of freedom. Such a treatment corresponds to nonequilibrium situations produced 
by external stirring mechanisms. (At finite temperatures, too, it forms the basis of the 
two-fluid model, wherein V pertains to the thermal cloud.) Thus, the center of mass of the 
condensate, in a local volume, moves relative to that of the interpenetrating cloud within 
with velocity V</> — V. Taking into account the local 'rigid-body' rotation associated with 
the cloud vorticity V x V, too, the kinetic energy density of the flow is, therefore, 



Here k is a parameter (with dimensions of moment of inertia per unit volume), which should 
depend on the macroscopic geometry of the flow. The above modifies the Gross-Pitaevskii 
energy functional (II]) via its last term, according to 




EW,il>*,V\= Jd 3 x |-(V v ^) 2 -^ o + |n o 2 + imn o (V0-V) 2 + ifc(Vx V) 2 
= jd 3 x [^i(V-^V)^| 2 -^| 2 +|M 4 +^(VxV) 2 . 
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This GL-like functional is invariant under local gauge transformations 
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The induced velocity V of the zero-point cloud, thus, plays the role of a gauge field with 
which the condensate particles interact, entailing a gauge interaction between the flows of 
the condensate and the cloud. In other words, the applied vorticity couples to the superflow 
in a gauge-field manner by means of inducing macroscopic internal vorticity V x V via 
the cloud. The disruption of the superflow by external vorticity is like the disruption of 
supercurrents by external magnetic fields. 

The field equations that follow from minimization of ([2]) are the well-known GL equations 

fcV x V x V= y(V>W*--0*VV>) -mV|<0| 2 (3) 

which replace the Gross-Pitaevskii equation of the standard condensate model. Their solu- 
tion subject to the flow conditions, in principle, determine ip and V for the flow. Note that 
in the zero-temperature limit we are considering, only the configuration of minimum energy 
contributes to the partition integral 

Z = J d[VMV>*MV] $[V]exp (-^m V]) (4) 

where the functional $[V] is a gauge fixing factor {5]. (It restricts the V-integrals to the 
physical values of the field V and thus renders these integrals finite. $ can be chosen from 
a variety of the gauge fixing factors (corresponding to the transverse (Landau) gauge, the 
Feynman gauge or the axial gauge for example), but the result will be independent of that 
choice, of course.) Hence, the GL equations flH]) can be considered exact here. 

It should be mentioned that the above GL functional has been already derived by Chelaflo- 

n n 

res who adopts the term "depletion" for the non-condensate fraction |8j. In his deriva- 
tion, he merely demands local gauge invariance for the Gross-Pitaevskii energy functional 
and identifies the emerging gauge field with the depletion velocity. However, he gives no 
physical justification. Moreover, the depletion velocity is not correlated with external agents 
and there is no mention of the source of this net macroscopic current (the microscopic mo- 
tions may be random after all, with no net current). Needless to say, without such physical 
ingredients the model cannot be applied to the study of superfluid turbulence. 
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IV. TRANSITION TO TURBULENCE 



The new gauge-field model of zero-temperature superfluidity, defined by the energy func- 
tional 02]), thus, coincides with the GL model of superconductivity in the superconductive 
phase. This explains the peculiar similarity, discussed above, between the two phenomena. 
Since \x > 0, the spontaneous breaking of the local gauge symmetry of ([2]) results in the 
Higgs mode. The Goldstone mode of the condensate model has transformed, giving way to 
a finite penetration depth for the induced internal velocity V and its vorticity field V x V. 
The exclusion of the interior vorticity, which occurs below a critical stirring force, is just the 
superfluid Meissner effect. The Meissner effect breaks down when the applied vorticity is 
too large. Superfluid flow can be divided into two classes according to how this breakdown 
occurs. This is determined by the GL parameter of the model k = y/gk/2h, which gives 
the relative strength of the short-range and gauge interaction terms, according to k < 1/ \/2 
(type I flow) or k > l/y/2 (type II flow). Thus, the flow types depend (through k) on the 
geometry of the flow, which is reconciling. In type I flows (such as the pipe flow), laminar 
flow is abruptly destroyed when the strength of the applied stirring force rises above a crit- 
ical value. In type II flows (such as the vessel flow), raising the applied vorticity above a 
first critical value leads to an intermediate mixed flow, in which an increasing amount of 
vorticity penetrates the superfluid through arrays of quantized vortices. At a second critical 
value, laminar flow is destroyed. 

To see how the above scenario for transition to turbulence arises, we apply the GL 
equations Q to a superfluid subject to external stirring. It is more convenient to start 
with a turbulent system and study the transition to laminar superflow flow instead. For 
sufficiently strong stirring, laminar superflow is destroyed. If we continuously decrease the 
applied vorticity, at a certain critical value, superflow begins to form. Two situations arise 
according to the two distinct solutions of the GL equations: 

(i) A condensate phase {ip ^ 0) appears abruptly in the bulk of the fluid, where is 
the same at all points. Working in the gauge ip =Real, the GL equations read: 

h 2 l 

V 2 ip + -mV 2 ip + gip 3 = flip 

2m 2 

fcV x V x V = -mVip 2 . (5) 
For ip = const., these imply the solution V = 0, ip — J fj,/g corresponding to a complete 
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expulsion of vorticity in the bulk. However, boundary effects are important here. Near the 
boundary, V and ip are no longer constant. Using ip = yjjjg as a zeroth approximation, 
equations (jSJ) yield V • V = and we obtain the London-like equation 

A 2 V 2 V = V 

where A = J~gkJJIm is the penetration length for the induced velocity/vorticity. The London 
equation explains the Meissner effect. The magnitude of the induced vorticity can be taken 
to be 20, where O is the applied angular velocity. This is valid because the superflow is weak 
near the boundary and the resulting corrections to the vorticity are, therefore, negligible. 
The induced cloud velocity near the boundary, then, has a magnitude V = OX. Since 
\mV 2 < \i (by equations K),we obtain a critical value, O c , for the applied angular velocity 
according to 

O < O c = S 2 ^-. 

As we decrease the applied angular velocity, Meissner effect takes place at O c . 

(ii) A condensate phase appears in the bulk of the fluid by spontaneous nucleation of 
condensate regions. In regions where nucleation occurs, superflow is just beginning to appear 
and, therefore, \tp\ is small. The GL equations can be linearized to give 

Also we can set V x V = 2fl. Here, this is valid because the current density of the condensate 
is of the order of l^j 2 and the resulting corrections to the vorticity are, therefore, negligible in 
the linear approximation. The above equation is, then, formally identical to the Shrodinger 
equation for a charged particle in a uniform magnetic field. Let us ignore boundary effects 
here, and consider an infinite domain. A condensate phase corresponding to an eigensolution 
of (jSJ) has an eigenvalue 

H = Knv 2 z + (2n + l)hO (7) 

where n is a nonnegative integer. It describes a particle with a constant velocity v z along f2 
and moving in a circle in the xy plane with angular frequency O. Therefore, the condensate 
phase which appears, in view of (j7|), for angular velocities O < 2 = fi/h, does not correspond 
to a complete expulsion of vorticity. 

In terms of the GL parameter, 2 = k\/20 c . When k < l/y/2 , we have ^2 < O c so that 
if we decrease O we first meet O c at which Meissner effect takes place. This corresponds 
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to type I flows. On the other hand when k > l/\/2, since fi 2 > the condensate phase 
does not correspond to Meissner expulsion for fl c < Q < f^. What we obtain here is that 
vorticity penetrates the bulk in the form of a dense array of quantized vortex lines, just as 
in superconductivity This corresponds to type II flows. 



V. THE WEAK-COUPLING LIMIT: LAMINAR FLOW 



For sufficiently weak stirring, the coupling of the induced internal vorticity to the su- 
perflow is weak and the flow is expected to be laminar. Corroboratively indeed, under this 
circumstance the two flows decouple and their gauge interaction reduces to an effective (long- 
range) Biot-Savart interaction among the condensate particles alone. This can be shown as 
follows. We rewrite the energy functional ([2]) as 

h 2 



E 



d x 



— iwr + {\mv 2 - »m 2 + fn 4 



fc(VxV) 2 -j .V 



(8) 



where 



ih 



Jo 



is the current density of the condensate. Under sufficiently weak stirring, such that \mV 2 <C 
/i, OH]) reduces to 



E = E + (Fx 



^(VxV) 2 -j 



V 



(9) 



where Eq is given by ([I]), of course. The gauge invariance is preserved because of the 
conservation of current density of the superflow (V • j = 0), which follows from the field 
equations derived from (Q. Since entropy is fixed at T « 0, the energy required by the 
condensate to exchange a particle with the cloud is given by its chemical potential //. The 
condition ^mV 2 fi, therefore, expresses that the induced kinetic energy of the cloud 
particles is not sufficient to disrupt the condensate by causing particle exchange, justifying 
the condensate current conservation. The two flows hence uncouple, which mathematically 
follows from substituting Q in (J3J) and integrating out the gauge field V from the partition 



function, according to [5( 

J d[V] $[V]exp 

oc exp 
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The result is 



J d[ip]d[ip* 



\exp[-—E[i/),if>* 



where 

E^,^] = EofaW -±-Uxj ^ jo(x)-j (xO (1Q) 

8nk J J |x — x'| 

is the effective energy. Again, the configuration of minimum energy only counts. Thus, 
under sufficiently weak stirring and as a by-effect of flow decoupling, there emerges a long- 
range Biot-Savart interaction between the condensate particles in addition to their short- 
range interaction. To minimize ([TO]) , the Biot-Savart term enforces local parallel ordering 
of currents, which is necessary to maintain laminar flow. Moreover, it may happen that 
vortex lines form, as mentioned before, to minimize E , in which case the Biot-Savart term 
favors their parallel ordering. This is possible only in appropriate geometries (corresponding 
to k > l/\/2) and explains the laminar vortex phase of type II flows. The Biot-Savart 
interaction between vortex lines is well-known in the context of the vortex filament model 
9|. Its automatic appearance here, as a by-effect of the induced gauge-field flow, is a 
corroborating feature of the model. 

To undergo transition to turbulent flow, the internal gauge flow must not decouple. The 
role of its induced vorticity is indispensable in destroying laminar flow, just as superconduc- 
tivity is destroyed by induced internal magnetization. 



VI. RECAPITULATION 



Due to strong interactions between superfluid particles, the fraction of particles in the 
zero-point cloud is large. When externally agitated, they contribute to the macroscopic 
flow via their induced average motion. Indeed, it is through the excitation of these internal 
motions that the applied vorticity couples to the superfluid. The internal relative motion 
of the flows of the cloud and the condensate defines a new degree of freedom, which is the 
essence of superfluidity in the zero-temperature limit. It entails a gauge coupling between 
the two interpenetrating flows and results in a GL-like model. The applied vorticity, thus, 
couples to the superflow in a gauge-field manner by means of inducing internal vorticity via 
the cloud. For sufficiently weak stirring, the applied vorticity cannot disrupt the condensate 
and the two flows decouple. The ensuing flow is a laminar flow of the condensate, where, as 
a by-effect of decoupling, a Biot-Savart interaction emerges between the condensate particles 
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in addition to their self-interaction. Here, two types of laminar flow may generally form: (i) 
The external vorticity is expelled from the fluid interior by the Meissner effect. The Biot- 
Savart interaction enforces local parallel ordering of currents and, thus, maintains laminar 
flow, (ii) The vorticity penetrates the fluid via formation of quantized vortices. The Biot- 
Savart interaction enforces parallel ordering of the vortex lines and, thus, forms vortex 
flow. The condition under which the second flow type is possible is determined by the GL 
parameter of the model, which, corroboratively, depends on the flow geometry. Finally, 
when the stirring is strong, the two flows are fully coupled and the laminar superflow is 
destroyed by the internal vorticity of the cloud. 
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